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1 Introduction 

High energy nucleus-nucleus collisions provide the means of creating nuclear 
matter in conditions of extreme temperature and density [H El Ej. The 
kinetic energy of the incident projectiles would be dissipated in the large 
volume of nuclear matter involved in the reaction. At large energy or baryon 
densities, a phase transition is expected from a state of nucleons containing 
confined quarks and gluons to a state of "deconfined" (from their individual 
nucleons) quarks and gluons, in chemical and thermal equilibrium, covering a 
volume that is many units of the confining length scale. This state of nuclear 
matter was originally given the name Quark Gluon Plasma (QGP) [I], a 
plasma being an ionized gas. 

A typical proposed phase diagram of nuclear matter [5j is shown in Fig. [I] 
together with the presumed trajectories of the evolution of the medium for 
collisions at RHIC and LHC cm. energies, where the axes are the temper- 
ature T vs. the baryon chemical potential \±b- The temperature for the 
transition from the Quark Gluon Plasma (QGP) to a hadron gas is taken as 
170 MeV for = and the phase boundary is predicted to be a smooth 
crossover down to a critical point below which the phase boundary becomes 
a first order phase transition. Also shown are idealized trajectories for the 
RHIC cm. energy scan and future experiments at FAIR which are being 
performed in order to find the QCD critical point. 
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Figure 1: (left) A proposed phase diagram for nuclear matter [5]: Temper- 
ature, T, vs Baryon Chemical Potential, \±b- 

2 Discovery of the QGP 

The QGP was discovered at RHIC, and announced on April 19, 2005. How- 
ever the results at RHIC [2] indicated that instead of behaving like a gas 
of free quarks and gluons, the matter created in heavy ion collisions at 
nucleon-nucleon cm. energy ^snn — 200 GeV appears to be more like a 
liquid. This matter interacts much more strongly than originally expected, 
as elaborated in peer reviewed articles by the 4 RHIC experiments [3 El El [9] , 
which inspired the theorists [10J to give it the new name "sQGP" (strongly 
interacting QGP). These properties were quite different from the "new state 
of matter" claimed in a press-conference [llj by the CERN fixed target 
heavy ion program on February 10, 2000, which was neither peer-reviewed 
nor published. 

In spite of not being published, the CERN press-release had a major ef- 
fect on the press in the United States resulting in an article on the front page 
of the New York Times [12]. Ironically, on this very same front page was an 
article announcing that the true version of the famous Italian sausage, Mor- 
tadella, would, for the first time, be allowed to be imported into the United 
States. A photograph of the iconic Bologna sausages appeared right next 
to the article about the CERN "qgp". Unfortunately, the first European 
Baloney to arrive in the U. S. was the CERN announcement. Q 

1 It is important for the reader of these proceedings to be aware that a high official 
of CERN was in the audience during this talk and made no objection to this comment. 
Furthermore, the author has a long, positive and productive relationship with this great 
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3 What's new at BNL and the RHIC machine this 
year. 

RHIC (Fig. § is one of only two Heavy Ion Colliders in the world, the other 
being the CERN LHC. RHIC is composed of two independent rings of su- 




Figure 2: Aerial view of RHIC/AGS facility. The two principal experiments 
still running are PHENIX and STAR. The LINAC is the injector for polar- 
ized protons. The TANDEM injector for Ions will be replaced by the Electon 
Beam Ion Source (EBIS) in fall 2011. The dome is the decommissioned High 
Flux Beam Reactor while the enclosed ring on the lower right center is the 
National Synchrotron Light Source (NSLS). 

perconducting magnets, 3.8 km in diameter (see Fig.[I^, below). RHIC can 
collide any species with any other species and so far has provided Au+Au, 
d+Au, Cu+Cu collisions at 12 different values of nucleon-nucleon cm. en- 
ergy y/sNN- For the past two years, an Au+Au energy scan has been per- 
formed with Energy/beam 3.85, 5.75, 9.8, 13.5, 19.5, 31.2, 100 GeV/nucleon. 
Also RHIC is the world's first and only polarized proton collider. The per- 
formance history of RHIC with A+A and polarized p-p collisions is shown in 
Fig. [3] At present RHIC operates at 15 times design luminosity for Au+Au 
and has shown a factor 2 progress in integrated luminosity C\ n t per week 
from Run-4 (2004) to Run-7 to Run-10 to Run-11. This year, 3 dimen- 
sional stochastic cooling was introduced for the Au+Au collisions which will 

laboratory and has praised its many successes. Thus he feels justified in commenting on 
one of their rare misjudgments. See reference [3 for a detailed scientific discussion. 
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Figure 3: a) (left) Au+Au performance, b) (right) polarized p-p perfor- 
mance. Courtesy Wolfram Fischer. 

improve the storage lifetime. Future improvements in the longitudinal pro- 
file, i.e. smaller diamond size without bucket migration, will be made with 
increased longitudinal focusing from new 56 MHz storage r.f. now under 
construction, with commissioning planned for Run-14. 

A significant improvement for the upcoming Run- 12 will be the replace- 
ment of the 40 year old Tandem Van de Graaff injector with an Electron 
Beam Ion Source (EBIS). A 10 A electron beam creates the desired charge 
state(s) in a trap within a 5 T superconducting solenoid. This is then ac- 
celerated through the RFQ and linac and injected into the AGS Booster 
(Fig. [5]). All ion species including noble gases, uranium and polarized 3 He 
are available. A uranium cathode has been received in preparation for a pilot 
U+U run in 2012. Commissioning of the EBIS started during early 2011 and 
it has already operated and supplied He + , He 2+ , Ne 5+ , Ne 8+ , Ar 11+ , Ti 18+ 
and Fe 20+ for the NASA Space Radiation Research Laboratory (NSRL) at 

bnl png. 

3.1 BNL's Superconducting Magnet Division 

All of the upgrades to the RHIC machine mentioned above involved super- 
conducting magnets built or developed in BNL's Superconducting Magnet 
Division. The Magnet Division also developed the RHIC machine super- 
conducting magnets (Palmer magnet [11]) which are the basis for the other 
post-Tevatron machines such as HERA and the LHC. The RHIC dipole 
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Figure 4: a) (left) Photo of the RHIC machine composed of two independent 
rings, with a total of 1740 superconducting dipole, quadrupole and corrector 
magnets, b) (right) Cross section of RHIC dipole. 

design (Fig. ^p) is based on a relatively large bore (80 mm inner diame- 
ter), single-layer "cosine theta" coil, wound from a (partially) keystoned, 
kapton-insulated, 30-strand Rutherford-type cable, arranged in coil blocks 
with intervening copper wedges, in order to meet the stringent field quality 
specifications, and mechanically supported by a laminated, cold steel yoke 
encased in a stainless steel shell. The shell contains the helium and is also 
a load bearing part of the assembly. 

It is amusing to note that in order to build the RHIC magnets for the 
purpose of making the QGP and studying its phase diagram, it is impor- 
tant to understand another phase-diagram, that of Fe+C, i.e. magnet steel 
(Fig. [5]). The Fe+C phase diagram is quite complicated with many phases, 
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Figure 5: Phase diagram of iron-carbon alloys [15J. 
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but is well known; while the proposed phase diagram of nuclear matter 
(Fig. [I]) seems much simpler, probably because it is largely unknown. 

3.2 Latest results from the Magnet Division 

BNL's superconducting magnet division is an international resource and ac- 
tively participates in many projects, several in the news this year. The 
CERN Courier of March 2011 featured on its cover (Fig. [6^1) the beautiful 
magnet used to trap anti- hydrogen for 1000 seconds in the ALPHA exper- 
iment [16] in the Antiproton Decelerator at CERN. The octupole magnet 




Figure 6: a) (left) CERN Courier cover March 2011. b) (right) Antihydrogen 
synthesis and trapping region of the ALPHA apparatus [16J. The atom-trap 
magnets, the modular annihilation detector and some of the Penning trap 
electrodes are shown (not to scale). 

(Fig. ^p) was built by the BNL Magnet Division and provides a very pure 
octupole field which is zero on the axis and rises sharply near the vacuum 
wall to keep the anti-hydrogen confined radially. 

Another recent press release featuring BNL magnets was the "Indication 
of Electron Neutrino Appearance" at the T2K (Tokai to Kamioka) Exper- 
iment in Japan on June 15, 2011 [18J. BNL provided 5 superconducting 
dipole corrector magnets in the proton beam at the Japan Proton Acceler- 
ator Research Complex (J-PARC) in Tokai which produced the /i-neutrinos 
that were detected in the Super-Kamiokande detector after transforming to 
e-neutrinos over the 295 km flight path. 

Other work of the Magnet Division includes: NbTi magnets for the RHIC 
Electron Lens upgrade; NbTi final focus quadrupole for the ILC; NbaSn 11.5 
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Tesla strand-test-barrel magnet and two coils for the LHC luminosity up- 
grade; High Temperature Superconductor (HTS) Quadrupole for the Facil- 
ity for Rare Isotope Beams at Michigan State University; Spare NbTi dipole 
for LHC; HTS solenoid R&D for muon collider and Energy Storage; Nb3Sn 
Open Midplane Dipole for the Muon Collider. 

3.3 National Synchrotron Light Source-II 

In addition to accelerators for high energy particle and nuclear physics, 
Brookhaven has been innovative in synchrotron radiation light sources. The 
National Synchrotron Light Source (NSLS) which started operations in 1982 
was the first to use the Chasman-Green double-bend achromat lattice [TT] , 
which is now the standard lattice at the major synchrotron light sources 
worldwide. A new third generation light source, NSLS-II, 4.66 times larger 
in circumference than NSLS, is now under construction at BNL, with unique 
design features of high brightness, small source size and long beam lines 
which will replace NSLS in 2014. NSLS-II is designed to deliver photons 
with high average spectral brightness in the 2 keV to 10 keV energy range ex- 
ceeding 10 21 ph/s/0.1%BW/ mm 2 /mrad 2 . The spectral flux density should 
exceed 10 15 ph/s/0.1%BW in all spectral ranges. This cutting-edge perfor- 
mance requires the storage ring to support a very high-current electron beam 
(I = 500 mA) with sub-nm-rad horizontal emittance (down to 0.5 nm-rad) 
and diffraction-limited vertical emittance at a wavelength of 1 A (vertical 
emittance of 8 pm-rad) [T9] . 

4 QGP Physics — Highlights from RHIC 

Given that I already said that the QGP was discovered at RHIC, what 
further studies are important? 

The QGP is the only place in the universe where we can in principle 
and in practice study Quantum Chromo-Dynamics (QCD) for color-charged 
quarks in a color-charged medium. For instance, how long will it take before 
we understand the passage of a quark through a QGP medium as well as we 
understand the passage of a muon through Copper in QED (Fig. [7])? 

Of course, in addition to understanding the behavior of QCD in a medium, 
the central goal of our field is a quantitative study of the phases of nuclear 
matter. This requires a broad, quantitative study of the fundamental prop- 
erties of the QGP including the extraction of the transport coefficients of 
the medium such as critical temperature, T c , speed of sound, c s , the ratio 
of shear viscosity to entropy density, 77/5, etc. To help understand how we 
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shall proceed to address these issues, it is important to understand how we 
got to this point. 

This is the 11th year of RHIC operation, and the two major detectors 
PHENIX and STAR which study the QGP at RHIC are basically first round 
detectors with a few incremental upgrades. Thus, the design of these de- 
tectors was heavily influenced by the c. 1990 view of the signatures of the 
QGP, which as noted above were quite different from what was discovered. 

4.1 J/^f suppression — the original "gold-plated" QGP signa- 
ture 

Since 1986, the 'gold-plated' signature of deconfinement was thought to 
be J/^f suppression. Matsui and Satz [21 j proposed that J/^f production 
in A+A collisions would be suppressed by Debye screening of the quark 
color charge in the QGP. The J/^f is produced when two gluons interact to 
produce a c, c pair which then resonates to form the J/^f. In the plasma 
the c, c interaction is screened so that the c, c go their separate ways and 
eventually pick up other quarks at the periphery to become open charm. 

"Anomalous suppression" of J/^f was found in ^snn — 17.2 GeV 
Pb+Pb collisions at the CERN SpS [22J (Fig. This is the CERN fixed 
target heavy ion program's main claim to fame: but the situation has al- 
ways been complicated because the J/\I/ is suppressed in p+A collisions. 
For example, in v / s]vjv = 38.8 GeV p+A collisions [23J (Fig. the Drell- 
Yan qq —> fi + fi~ cross-section per nucleon is constant as a function of mass 
number, A, which indicates the expected absence of shadowing in a nucleus 
for point-like production processes; while the J/\& and T cross sections per 
nucleon are suppressed by an amount A a with a = 0.920 ± 0.008 for both 
J/\& and and a = 0.96 + 0.01 for both the T ls and T 2s+3s . This is called 
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Figure 8: a) (left) Total cross section for J /iff production divided by AB in A+B 
collisions at 158-200 A GeV [22] (note the *). b) (right) A dependence of charmo- 
nium and Drell-Yan pair production in 800 GeV p+A collisions [23] expressed as 
the ratio of heavy nucleus to deuterium cross sections per nucleon. The dashed 
lines are fits to A a for the CNM effect, with the values of a indicated. 



a Cold Nuclear Matter or CNM effect and is shown as the line with a = 0.92 
on Fig. [8^i. The "Anomalous suppression" is the difference between the data 
point at AB = 208 2 and the line, provided that the CNM effect is the same 
at y/^T = 17.2 and 38.8 GeV. 

The search for J /iff suppression and thermal photon/dilepton radiation 
from the QGP drove the design of the RHIC experiments. 



4.2 Detector issues in A+A compared to p-p collisions 

A main concern of experimental design in RHI collisions is the huge mul- 
tiplicity in A+A central collisions compared to p-p collisions. A schematic 
drawing of a collision of two relativistic Au nuclei is shown in Fig. [9^i. In 
the center of mass system of the nucleus-nucleus collision, the two Lorentz- 
contracted nuclei of radius R approach each other with impact parameter 
b. In the region of overlap, the "participating" nucleons interact with each 
other, while in the non-overlap region, the "spectator" nucleons simply con- 
tinue on their original trajectories and can be measured in Zero Degree 
Calorimeters (ZDC), so that the number of participants can be determined. 
The degree of overlap is called the centrality of the collision, with 6^0, 
being the most central and b ~ 2i?, the most peripheral. The maximum 
time of overlap is r Q = 2R/^c where 7 is the Lorentz factor and c is the 
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Figure 9: a) (left) Schematic of collision of two nuclei with radius R and 
impact parameter b. The curve with the ordinate labeled da/dn^ represents 
the relative probability of charged particle multiplicity n c h which is directly 
proportional to the number of participating nucleons, Af part . b) (right) Trans- 
verse energy (Et) distribution in Au+Au and p-p collisions at y / s]v7v = 200 
GeV from PHENIX [25]. 



speed of light in vacuum. 

The energy of the inelastic collision is predominantly dissipated by mul- 
tiple particle production, where iV c h, the number of charged particles pro- 
duced, is directly proportional [9J to the number of participating nucleons 
(iVp ar t) as sketched on Fig.[9^i. Thus, iV c h or the total transverse energy Ej 1 
in central Au+Au collisions is roughly A times larger than in a p-p collision, 
as shown in the measured transverse energy spectrum in the PHENIX de- 
tector for Au+Au compared to p-p (Fig. [9b) and in actual events from the 



STAR and PHENIX detectors at RHIC in Fig. 10 The impact parameter 




Figure 10: a) (left) A p-p collision in the STAR detector viewed along the 
collision axis; b) (center) Au+Au central collision at ^snn — 200 GeV in 
STAR; c) (right) Au+Au central collision at y/s^ = 200 GeV in PHENIX. 
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b can not be measured directly, so the centrality of a collision is defined in 
terms of the upper percentile of n c h or Et distributions, e.g. top 10%-ile, 
upper 10 — 20%-ile. Unfortunately the "upper" and "-ile" are usually not 
mentioned which sometimes confuses the uninitiated. 

In Fig. [TTJ measurements of the charged particle multiplicity density 
dNch/dr] at mid-rapidity, \rj\ < 0.5, relative to the number of participating 
nucleons, iV par t, are shown as a function of centrality for ^s^n = 200 GeV 
Au+Au collisions at RHIC [25J together with new results this year from 
ALICE in = 2.76 TeV Pb+Pb collisions at LHC [28J. The results are 

expressed as (dN&i/drj) / (Af part /2) for easy comparison to p-p collisions. 
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Figure 11: Dependence of (dN^/drf) /(iV par t/2) on the average number of 
participants (Af part ) in bins of centrality, for Pb+Pb collisions at ^/snn = 
2.76 TeV [28] and Au+Au collisions at ^/s^ = 0.200 TeV [25]. The scale 
for the lower-energy data (right side) differs by a factor of 2.1 from the scale 
for the higher-energy data (left side). 



The LHC data show the effect well known from RHIC that dN^/drj does 
not depend linearly on A^ part , since {dN^/drj) /(^V part /2) is not a constant 
for all Af part . However the data also show the amazing effect that the ratio 
of (dNch/ drj) / (N part / 2) from LHC to RHIC is simply a factor of 2.1 in every 
centrality bin. Thus the LHC and RHIC data lie one on top of each other 
by simple scaling of the RHIC measurements by a factor of 2.1. This is an 
incredibly beautiful result which shows that in going from p-p to A+A col- 
lisions, the charged particle production is totally dominated by the nuclear 
geometry of the A+A collisions represented by the number of participating 
nucleons, A^ part , independently of the nucleon-nucleon cm. energy, ^/snn- 

Since it is a huge task to reconstruct the momenta and identity of all the 
particles produced in these events, the initial detectors at RHIC [26J con- 
centrated on the measurement of single-particle or multi-particle inclusive 
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variables to analyze RHI collisions, with inspiration from the CERN ISR 
which emphasized those techniques before the era of jet reconstruction [27] . 
There are two major detectors in operation at RHIC, STAR and PHENIX, 
and there were also two smaller detectors, BRAHMS and PHOBOS, which 



have completed their program. As may be surmised from Fig. 10, STAR, 
which emphasizes hadron physics, is most like a conventional general pur- 
pose collider detector, a TPC to detect all charged particles over the full 
azimuth (Ac/) = 2tv) and ±1 units of pseudo-rapidity (77); while PHENIX 
is a very high granularity high resolution special purpose detector cover- 
ing a smaller solid angle at mid-rapidity, together with a muon-detector at 
forward rapidity [29] . 

One nice feature of the STAR detector is the ability to measure the 
mass/charge of a particle from its momentum/charge and time of flight, 
and then use dE/dx measured in the TPC to determine the charge. In this 
way STAR has observed many anti-nuclei, notably this year the "Observa- 
tion of the antimatter helium-4 nucleus" in Au+Au collisions [30] (Fig.[l2^i). 
The differential invariant yields d 2 N /(2tv pTdpTdy) per central Au+Au col- 
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Figure 12: a) (left) Number of standard deviations, n adE/dx , of dE/dx reso- 
lution from the expected value for 4 He, for negative and positive particles as 
a function of calculated mass, and projected counts for 
b) Differential invariant yields as a function of baryon number, B. 



"2 < n <r d E/ dx < 3 - 



lision at -^/snn = 200 GeV as a function of baryon number B, evaluated at 



Pt/\B\ = 0.875 GeV/c, are shown in Fig. 12 3 and show a steady exponential 
decrease with increasing B [31] . The anti- nuclei are made by coalescence of 
the large number of n and p produced, an advantage of the high multiplicity. 

PHENIX is designed to measure and trigger on rare processes involving 
leptons, photons and identified hadrons at the highest luminosities with the 
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special features: i) a minimum of material (0.4% X Q ) in the aperture to avoid 
photon conversions; ii) possibility of zero magnetic field on axis to prevent 
de-correlation of e + e _ pairs from photon conversions; iii) Electro-Magnetic 
Calorimeter (EMCal) and Ring Imaging Cherenkov Counter (RICH) for 

identification and level- 1 trigger; iv) a finely segmented EMCal (£77, 
5<j) — 0.01 x 0.01) to avoid overlapping showers due to the high multiplicity 
and for separation of single-7 and 7r° up to pt ~ 25 GeV/c; v) EMCal 
and precision Time of Flight measurement for particle identification. Some 
results uniquely possible with this detector such as measurements of direct 
photons via internal conversion to e + e~ pairs will be discussed below. 

In addition to the large multiplicity, there are two other issues in RHI 
physics which are different from p-p physics: i) space-time issues, both in 
momentum space and coordinate space — for instance what is the spatial 
extent of fragmentation? is there a formation time/distance?; ii) huge az- 
imuthal anisotropics of particle production in non-central collisions (collo- 
quially collective flow) which are very interesting in their own right and 
provide much richer features than originally envisaged. 



5 Collective Flow 



A distinguishing feature of A+A collisions compared to either p-p or p+A 
collisions is the collective flow observed. This effect is seen over the full range 
of energies studied in heavy ion collisions, from incident kinetic energy of 
100 A MeV to cm. energy of ^/s^ = 200 GeV [32J. Collective flow, or 
simply flow, is a collective effect which can not be obtained from a superpo- 
sition of independent N-N collisions. Immediately after an A+A collision, 
the overlap region defined by the nuclear geometry is almond shaped (see 



Fig 13) with the shortest axis along the impact parameter vector. Due to the 



reaction plane breaking the (f) symmetry of the problem, the semi-inclusive 
single particle spectrum is modified by an expansion in harmonics |35j of the 
azimuthal angle of the particle with respect to the reaction plane, (j)—^R [36] , 
where the angle of the reaction plane is defined to be along the impact 
parameter vector, the x axis in Fig. [13) 



Ed 3 N 
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(1) 



The expansion parameter ^2, called elliptical flow, is predominant at mid- 
rapidity. In general, the fact that flow is observed in final state hadrons 
shows that thermalization is rapid, so that hydrodynamics comes into play 
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Figure 13: (left) Almond shaped overlap zone generated just after an A+A 
collision where the incident nuclei are moving along the ±z axis. The reac- 
tion plane by definition contains the impact parameter vector (along the x 
axis) [33J. (right) Measurements of elliptical- flow (i^) for identified hadrons 
plotted as divided by the number of constituent quarks n q in the hadron 
as a function of (a) pr/^qi (b) KET/n q [3i] . 



at a time, To, which is before the spatial anisotropy of the overlap almond 
dissipates. At this early stage hadrons have not formed and it has been 
proposed that the constituent quarks flow [37J, so that the flow should be 
proportional to the number of constituent quarks n q: in which case V2/n q as a 
function oipT/n q would represent the constituent quark flow as a function of 
constituent quark transverse momentum and would be universal. However, 
in relativistic hydrodynamics, at mid-rapidity, the transverse kinetic energy, 
tut — mo = (tt — 1)ttio = KEt, rather than pt, is the relevant variable; 
and in fact V2/n q as a function of KET/n q seems to exhibit nearly perfect 
scaling [34] (Fig.[l3]b). 

The fact that the flow persists for pt > 1 GeV/c (Fig. [l4^ i) implies that 
the viscosity is small [38], perhaps as small as a quantum viscosity bound 
from string theory [39], rj/s = 1/(Att) where 77 is the shear viscosity and s 
the entropy density per unit volume. This has led to the description of the 
"sQGP" produced at RHIC as "the perfect fluid" [10J. An estimate [40] 
of r]/s for nuclear matter and for several common fluids, as a function of 
the fractional difference of the temperature from the critical temperature, 
at fixed pressure, is shown in Fig. 14 j. This particular estimate [ID] for 
the QGP at RHIC is quite close to the quantum bound (solid line). Also, 
empirically, for all common fluids rj/s is a minimum at or near the critical 
point [H] which might suggest that the conditions at RHIC energies are 
near the QCD critical point. 
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Figure 14: a) (left) Teaney's [38] predictions for V2(pt) for ideal (T s /to = 
0) and viscous hydrodynamics, where T s = |^ is the sound attenuation 
length, and To is the thermalization time, b) rj/s for various fluids at fixed 
pressure as a function of temperature T, where Tq is the temperature of the 
critical point of the liquid-gas phase transition [^Of ET]. 



5.1 Two-Particle Correlations and Flow 

In addition to measuring flow by the correlation of individual particles to 
the reaction plane, it is also possible to measure flow by the correlation of 
two-particles to each other. The advantage of this method is that one does 
not have to determine the reaction plane. Thus if two particles A and B are 
correlated to the reaction plane, but not otherwise correlated to each other, 
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then the correlation to the reaction plane induces a correlation of these two 
particles to each other which can be measured without knowledge of the 
reaction plane, 
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In p-p collisions there is no collective flow but there are strong two- 
particle azimuthal correlations due to di-jet production in hard-scattering 
(see Fig. [TO^t), which also exist in A+A collisions but are obscured by the 
large multiplicity [e.g. can you find a jet in Fig. [l0fc>] . Before the discovery 
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of jets, two-particle correlations were used extensively at the CERN-ISR to 
study hard-scattering (the production of particles with large transverse mo- 
mentum) in p-p collisions which was discovered there [27]. Due to the huge 
multiplicities in Au+Au collisions at RHIC, where for central Au+Au colli- 
sions at y/sNN = 200 GeV, there is an estimated nAr 2 x ^^f^ 1 ~ 375 GeV 
of energy in one unit of the nominal jet-finding cone, Ar = \J (At/) 2 + (A0) 2 , 
two-particle correlations were used exclusively for the first 10 years to study 
hard scattering at RHIC. 

Typical examples of the di-hadron measurements in p-p and Au+Au 
central (0-20%) collisions at ^/s^ = 200 GeV are shown in Fig. [Ts] [421 S3] 
which are presented as azimuthal distributions of the conditional yields of 
associated particles, with pT a , with respect to trigger particles with 3 < 
PT t < 10 GeV/c. The di-jet structure in p-p collisions is clearly indicated 




Figure 15: a-h) (left) PHENIX [32] azimuthal correlation conditional yield of 
associated h ± particles with pT a for trigger h ± with pT t for the various pT t ® 
PT a combinations shown, i) (right)- (top) PHENIX [43J azimuthal correlation 
function C(A0) of h ± with 1 < pT a < 2.5 GeV/c with respect to a trigger h ± 
with 2.5 < pT t < 4 GeV/c in Au+Au central collisions, where the line with 
data points indicates C(A(j>) before correction for the azimuthally modulated 
(^2) background, and the other line is the V2 correction which is subtracted 
to give the jet correlation function J(A0) (data points), j) (right)- (bottom) 
PHENIX D parameters [33] > the angular distance of the apparently displaced 
peak of the J(A0) distribution from the angle A0 = tv as a function of 
centrality, represented as the number of participants Af part , for the systems 
and cm. energies indicated. 
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by the gaussian-like strong azimuthal correlation peaks on the same side 
(Ac/) = <j) a — (j>t ~ 0) and away side {A(j> ~ tv rad.) relative to the trigger 
particle for all ranges of pT t and pT a measured. However, one of the many 
interesting features in Au+Au collisions is that the away side azimuthal jet- 
like correlations (Fig. [l5fc ) are much wider than in p-p collisions and show a 
two-lobed structure ("the shoulder" (SR)) at lower px t with a dip at 180°, 
reverting to the more conventional structure of a peak at 180° ("the head" 
(HR)) for larger pT t - 

The wide away-side correlation in central Au+Au collisions is signifi- 
cantly obscured by the large multiparticle background which is modulated 
in azimuth by the V2 collective flow of a comparable width to the jet corre- 
lation (Fig. [T5j) . After the V2 correction, the double peak structure ~ ±1 
radian from 7r, with a dip at tt radians, becomes evident. The double-peak 
structure may indicate a reaction of the medium to a passing parton in anal- 
ogy to a "sonic boom" or the wake of a boat, which was given the name 
"Mach Cone" [33], and has been under active study both theoretically [35] 
and experimentally. PHENIX characterizes this effect by the half-width D 
(~ 1.1 radian) of the Jet function, J(A0), the angular distance of the dis- 
placed peak of the distribution from the angle Acj) = tt. One of the striking 
features of the wide away side correlation is that the width D (Fig. [l5] j) 
does not depend on centrality, angle to the reaction plane, pT a and ^/sNN, 
which always seemed problematic to me if the effect were due to a reac- 
tion to the medium. Another suspicious issue is that the same effect occurs 
even for auto-correlations of particles with very low pt between 0.2 and 0.4 
GeV/c where any effect of hard-scattered partons should be submerged by 
the predominant soft physics (Fig. [16^) [35] . 

In addition to the Head/Shoulder or Mach Cone effect in two-particle 
correlations on the away-side, same-side correlations also show a new effect 
in A+A collisions called "the Ridge" [37]. This is seen in two-dimensional 
correlations in A77, Acf) (Fig. (16)3) [351 139] where the associated yield dis- 
tribution can be decomposed into a narrow jet-like peak at small angular 
separation which has a similar shape to that found in p-p collisions, and 
a component that is narrow in Acj) but depends only weakly on A 77, the 
"ridge." However, new results this past year have dramatically changed this 
picture. 

5.2 Triangular flow, odd harmonics 

For the first 10 years of RHIC running, and dating back to the Bevalac, all 
the experts thought that the odd harmonics in Eq. [T] would vanish by the 
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Figure 16: a) (left) Low pt like-sign pair azimuthal correlation function for 
0-5% central Au+Au collisions at ^snn = 200 GeV from charged hadrons 
with 0.2 < p Tl ,PT 2 < 0.4 GeV/c [46J. b) (right) "The Ridge" [M\M- 



symmetry (p —> (p + tt of the almond shaped overlap region [50] (Fig. 13). 
However, last year, 2010, an MIT graduate student an his Professor in ex- 
perimental physics, seeking (at least since 2006) how to measure the fluctu- 
ations of V2 in the PHOBOS experiment at RHIC, realized that fluctuations 
in the collision geometry on an event-by-event basis, i.e. the distribution 
of participants from event-to-event, did not respect the average symmetry. 
This resulted in what they called "participant triangularity" and "triangular 
flow" , or v% in Eq. [TJ which they measured using both PHOBOS and STAR 
data [51J0 

Many experiments presented measurements of v% at Quark Matter 2011 



this year, e.g. Fig. 17 [53J, and it was one of the most exciting results of the 



past year. There are two striking observations from Fig. 17 which indicate 
that fluctuations of the initial collision geometry are driving the observed 
^3: i) the centrality dependence of vsipr) is weak as one would expect 
from fluctuations, but i^Opt) which is most sensitive to the geometry of 
the "almond" -shaped overlap region tracks the change in eccentricity with 
centrality; ii) for the most central collisions (0-10%), where the overlap 
region is nearly circular so that all the v n are driven by fluctuations, i^Cpt)? 
v 3(pt), ^a{pt) are comparable. The fact that the observed collective flow 
of final state particles follows the fluctuations in the initial state geometry 



2 It was pointed out by Leticia Palhares in the discussion that a Brazilian group showed 
in 2009 that that the ridge and the cone, i.e. V3, does appear in an event- by-event 
hydrodynamics calculation without jets [52], but the MIT group [51] was the first to show 
it with real data. 



18 



0-10% 


10-20% 


20-30 % 


30-40 % 


40-50 % 


50-60 % 




Au+Au 200GeV 

• v 2{^ 2 } 
A V 3^ 3 } 
□ v 4 {^ 4 } 


it 


f 


if 
if 

f 


• 

• 




•' 

t 

i 




W 






/ y XL- 




i„„i„„i„„i„„i„„i„„i„„i„ 













0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 3 3.5 0.5 1 1.5 2 2.5 33.5 

P T [GeV/c] 



Figure 17: PHENIX [53] measurements of the v n parameters using Eq. [I] 
(with the appropriate reaction plane) as a function of pt for different cen- 
trality slices in ^/snn = 200 GeV Au+Au collisions. 

points to real hydrodynamic flow of a nearly perfect fluid (and convinces this 
author of the validity of hydrodynamics in RHI collisions, of which he was 
quite skeptical two years ago [27]). It is evident that ^3, a cos3(A0) term 
with lobes at = 0, 2tv/3 and 47r/3 « 0,2,4 radians, would explain the 
peaks at tt ± D radian in the two-particle correlations (Fig. [l5|j) and also 
why D^l radian independent of centrality and kinematic variables; while 
the lobe at = explains the ridge (Fig. (16)3). There is presently lots of 
activity to confirm in detail whether taking account of the odd harmonics in 
addition to V2 and V4 in the background of Fig. [15} will result in gaussian- 
like away-jet peaks in A+A collisions and the disappearance of the same-side 
ridge. 

6 RHIC beam energy scan — In search of the crit- 
ical point 

The past year has seen runs at RHIC with Au+Au collisions at cm. energies 
y/sjsrN — 7.7, 11.5 and 39 GeV, in addition to previous runs at 19.6, 62.4, 130 
and 200 GeV, to search for the onset of large fluctuations which should occur 
near a critical point in the phase diagram (Fig. [TJ. Such fluctuations in the 
K/tt and K/p ratio had been claimed at the CERN SPS fixed target program 
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near ^snn = 8 GeV and were presented [51] as "evidence of the onset of the 
deconflnement phase transition". At QM2011, STAR [55] presented many 



excellent results on this subject of which I show a small selection in Fig. 18 
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Figure 18: a) (left) iV c h distribution in the STAR detector for 5 values of 
■\/snn [55] : b) (center) K/tt ratio vs y/sjsnsr [55]; c) (right) Event-by-event 
fluctuation of K/p ratio [55] . 



Fig. 18 a, shows the multiplicity distribution dN^/dr] which maintains the 
characteristic Nuclear Geometry dominated shape (as in Fig. |9|, stretching 
to higher multiplicity, dN^/drf, with increasing ^/snn. Fig. Il8b shows the 
K + /tv + and K~ /tv~ ratios over the entire range of ^/s^n measurements. 
The maximum of the K + /tt + ratio near ^sjyjy — 8 GeV is naturally ex- 
plained [56] by the change in dominant K + production from K + A to K + K~ 
whose smooth increase with ^/snn can be seen from the K~ /tv~ ratio. The 
famous "horn" , or apparent discontinuity, at ^/snn = 8 GeV from the SPS 



data [54J is greatly smoothed when the new STAR data are added. Fig. 18 



shows a smooth variation of the the STAR measurements of the fluctua- 
tions of the event-by-event K/p ratio as a function of ^/saw, which differs 
dramatically from the huge effect claimed by the SPS Fixed Target mea- 
surements below 12 GeV, notably the change from negative to positive [5^] . 
There is no doubt in this author's mind that one must prefer the collider 
measurements, where the detector position at mid-rapidity in the cm. sys- 
tem is constant for all values of y/sjsrN, to the fixed target measurements, 
where the rapidity of the cm. system moves dramatically with respect to 
the the detector as ^snn varies. This is a major strength of the RHIC 
Beam Energy Scan program. 
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6.1 A press release during the school 

On June 23, 2011, shortly before I was to give these lectures, a press re- 
lease from LBL arrived claiming that "By comparing theory with data from 
STAR, Berkeley Lab scientists and their colleagues map phase changes in 
the QGP" [57]. Since I had criticized "physics by press-release" concerning 
the discovery of the QGP (above), I felt that I was obliged to review the 
physics behind this latest example, presumably a "Highlight from RHIC" . 

The subject is "Fluctuations of conserved quantities", in this case the 
net baryon distribution taken as p — p. Since there can be no fluctuations 
of conserved quantities such as net charge or net baryon number in the full 
phase space, one has go to small intervals [59] to detect a small fraction of 
the protons and anti-protons which then fluctuates, i.e. varies from event to 
event. The argument is that, e.g. the fluctuation of one charged particle in 
or out of the considered interval produces a larger mean square fluctuation of 
the net electric charge if the system is in the hadron gas phase with integral 
charges than for the QGP phase with fractional charges. 

However, while there are excellent statistical mechanical arguments about 
the utility of fluctuations of conserved quantities such as net baryon number 
as a probe of a critical point [58J, there are, so far, no adequate treatments of 
the mathematical statistics of the experimental measurements. Theoretical 
analyses tend to be made in terms of a Taylor expansion of the free energy 
F = —ThiZ around the critical temperature T c where Z is the partition 
function, or sum over states, which is of the form Z oc e~ ( " E ~^i^ i( ^ i ^ kT and 
\ii are chemical potentials associated with conserved charges Qi [58]. The 
terms of the Taylor expansion, which are obtained by differentiation, are 
called susceptibilities, denoted x- The only connection of this method to 
mathematical statistics is that the cumulant generating function in mathe- 
matical statistics is also a Taylor expansion of the In of an exponential: 

(3) 

t=0 

Thus, the susceptibilities are cumulants in mathematical statistics terms, 
where, in general, the cumulant K m represents the m th central moment with 
all m-fold combinations of the lower order moments subtracted^ For in- 
stance, K2 {(x — n) 2 ) = <T 2 , ft3 ((x — /x) 3 ), ^4 ((x — n) A ) —3^5 

^5 = ((x — /x) 5 ) — 10^3^2, where n — (x). Two so-called normalized or 

3 Note that factorial cumulants, also known as "Mueller moments" 60 , have been used 
previously for multiplicity distributions. 
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standardized cumulants are common in this field, the skewness, S = 
and the kurtosis, n — = ((x — fi) 4 ) /a 4 — 3. 

A sample |61j of STAR measurements of the distribution of net-protons 
in Au+Au collisions in the small interval 0.4 < pt < 0.8 GeV/c, \y\ < 0.5 



for different s/snn is shown in Fig. 19 a,. 
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Figure 19: a) (top-left) STAR [61] distribution of event-by-event p — p at 3 
values of ^/snn; b) (top-right) STAR published [62] measurements of ftcr 2 ; 
c) (bottom-left) Measurements from (b) as shown in Ref. [63] compared 
to the predicted ratio of susceptibilities (open crosses); d) (bottom-right) 
compilation [61 J of STAR measurements of no 1 . 

The moments in the form no 1 — ^4/^2 are shown from a previous STAR 



publication [62] in Fig. 193 while a plot, alleged to be of this same data, 
presented in the Lattice QCD theory publication that generated the press- 

[63] ; and a plot of the no 2 from the data of 
combined with the results from Fig. 19 3, is shown in Fig. [l9|d [61J. 



release, is shown in Fig. [19 
Fig 
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There are many interesting issues to be gleaned from Fig. [19 



The data point at 20 GeV in Fig. 19: is not the published one from (b), 
as stated in the caption [63], but the one from (d), which is different and 
with a much larger error. This, in my opinion, makes the data point look 



22 



better compared to the predicted discontinuous value of kg 2 — — 1.5 for the 
critical point at 20 GeV (open crosses) in contrast to the predictions of 1.0 
for both 62.4 and 200 GeV. The published measurements in (b) together with 
the new measurements in (d) are all consistent with kg 2 — 1; but clearly 
indicate the need for a better measurement at ^/snn = 20 GeV. Apart 
from these issues, the main problem of comparing Lattice QCD "data" to 
experimental measurements is that it is like comparing peaches to a fish, 
since the prediction is the result of derivatives of the log of the calculated 
partition function of an idealized system, which may have little bearing 
on what is measured using finite sized nuclei in an experiment with severe 
kinematic cuts. Maybe this is too harsh a judgement; but since this is the 
first such comparison (hence the press release), perhaps the situation will 
improve in the future. 

When I first saw the measured distributions in Fig. 19 a,, my immediate 
reaction was that STAR should fit them to Negative Binomial distributions 
so that they would know all the cumulants. However, I later realized that 
my favorite 3 distributions for integer random variables, namely, Poisson, 
Binomial, and Negative Binomial, are all defined only for positive integers, 
while the number of net-protons on an event can be negative as well as 
positive, especially at higher cm. energies. This is why somebody should 
work out the mathematical statistics of the net proton distribution as we 
did in PHENIX for the distribution of the difference in foreground (op- 
posite charge) and background (like charge) di-lepton events when both 
are Poisson distributed [64J. Until then, it is instructive to compare the 
values of kg 2 in Fig. [l9] to those from the well-known distributions: Pois- 
son, kg 2 — 1; Binomial, kg 2 — 1 — 6p + 6p 2 < 1; Negative Binomial, 
na 2 = 1 + 6/i/fc + 6/i 2 //c 2 > 1; Gaussian, all cumulants=0 for k > 2, so 
kg 2 — 0. The data favor Poisson (i.e. no correlation) everywhere, with 
some hint of Binomial. Nevertheless, if a future measurement would show 
a significant huge discontinuity of kg 2 similar to the theoretical prediction 
at -^/snn = 20 GeV, then even I would admit that such a discovery would 
deserve a press release, maybe more! 



7 Hard Scattering as a probe of the QGP 

One of the best probes found at RHIC to study the QGP is the hard- 
scattering of quarks and gluons (the constituents of the nucleon) which can 
be observed via inclusive single particle production at large transverse mo- 
mentum (pt) or by two-particle correlations with a high pt trigger. The 
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hard-scattering takes place in the initial collision of the highly Lorentz con- 
tracted nuclei. The scattered partons which emerge near 90° to the collision 
axis (the sweet-spot for such observations) pass through the medium formed 
and then fragment to jets of particles which are detected (Fig. 20). The most 
likely 

Y constituent reactions are: g+g -+ g+g, g + q — > g+q, 

q+q — » g + g, and g + q —> 7 + g (shown). This last re- 
action [55] . called direct-7 production or the inverse 
QCD Compton effect, is the most beautiful because 
the 7-ray participates directly in the hard scattering, 
then emerges from the medium without interacting 
and can be measured with high precision. No frag- 
mentation is involved and, in fact, these direct-7 are 
distinguished from e.g. 7 rays from fragmenting par- 
tons because they are isolated, with no accompany- 
ing fragments. Triggering on a direct-7 of a given pr t 
provides a 'beam' of partons (8/1 u-quark) with ex- 
actly opposite (thus well-known) initial pt = —PT t , 
so that the effect of the medium can be determined by measuring, for in- 
stance, the ratio of the transverse momentum pT a of the jet from the away- 
parton to that of the direct-7 trigger, denoted Xh = PT a /PT t , or equivalently, 
the fractional jet imbalance, 1 — x^, as used by CMS at LHC [661. 

Since hard-scattering at high pt > 2 GeV/c is point-like, with distance 
scale 1/pt < 0.1 fm, the cross section in p+A (A+A) collisions, compared to 
p-p, should be larger by the relative number of possible point-like encounters, 
a factor of A (A 2 ) for p+A (A+A) minimum bias collisions. When the 
impact parameter or centrality of the collision is defined, the proportionality 
factor becomes (Taa), the average overlap integral of the nuclear thickness 
functions. 



Figure 20: g + q —> 
7 + q reaction with 
outgoing 7 and q 
passing through the 
medium 



7.1 Jet quenching — suppression of high p T particles 

The discovery, at RHIC [57], that 7r°'s produced at large transverse mo- 
menta are suppressed in central Au+Au collisions by roughly a factor of 
5 compared to point-like scaling from p-p collisions is arguably the ma- 
jor discovery in Relativistic Heavy Ion Physics. For tt° (Fig. 21 1) [68] the 
hard-scattering in p-p collisions is indicated by the power law behavior p^ n 
for the invariant cross section, Ed 3 a /dp 3 , with n = 8.1 ± 0.05 for pt > 3 
GeV/c. The Au+Au data can be characterized either as shifted lower in 
energy relative to the point-like scaled p-p data, or down in magnitude, 
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Figure 21: a) (left) Invariant yield of tt° at y/^T = 200 GeV as a function 
of transverse momentum pt in p-p collisions multiplied by (Taa) f° r Au+Au 
central (0-10%) collisions compared to the Au+Au measurement [68J. b) 
(right) Raa(pt) for all identified particles so far measured by PHENIX in 
Au+Au central collisions at ^snn = 200 GeV. 



i.e. suppressed. In Fig. |21p , the suppression of the many identified parti- 
cles measured by PHENIX at RHIC is presented as the Nuclear Modification 
Factor, Raa(pt), the ratio of the yield of per central Au+Au collision (upper 
10%-ile of observed multiplicity) to the point-like-scaled p-p cross section: 

d 2 Nl A /dp T dyN A A m 
Raa(pt)= {TAA)d ^ dpTdy ■ (4) 

The striking differences of Raa(pt) in central Au+Au collisions for the 
many particles measured by PHENIX (Fig. (^lfc) illustrates the importance 
of particle identification for understanding the physics of the medium pro- 
duced at RHIC. Most notable are the equal suppression by a constant factor 
of 5 of 7T° and 77 for 4 < pt < 15 GeV/c, the equality of suppression of 
direct-single e ± (from heavy quark (c, b) decay) and 7r° at Pt^5 GeV/c, 
the non-suppression of direct-7 for pt > 4 GeV/c and the exponential rise 
of Raa of direct-7 for pt < 2 GeV/c [69J, which is totally and dramat- 
ically different from all other particles and attributed to thermal photon 
production by many authors (e.g. see citations in reference [69]). For pr^4 
GeV/c, the hard-scattering region, the fact that all hadrons are suppressed 
but direct-7 are not suppressed, indicates that suppression is a medium ef- 
fect on outgoing color-charged partons. The suppression can be explained by 
the energy loss of the outgoing partons in the dense color-charged medium 
due to coherent Landau-Pomeranchuk-Migdal radiation of gluons, predicted 
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in QCD [70J, which is sensitive to properties of the medium. Measurements 
of two-particle correlations (discussed below) confirm the loss of energy of 
the away-jet relative to the trigger jet in Au+Au central collisions com- 
pared to p-p collisions. However, as we shall see, lots of details remain to 
be understood. 

The most interesting new results this year concern: i) direct-7-hadron 
correlations at large pt to measure the fragmentation function in p-p colli- 
sions and to search for a possible modification of the fragmentation function 
in Au+Au; ii) 7r°-hadron correlations to compare the fractional jet imbal- 
ance at RHIC to the new LHC measurement [66] ; iii) measurement of the 
flow (^2) of direct-7 at low and high pt] and iv) confirmation at the LHC 
of the suppression of heavy quarks comparable to that of light quarks. 

7.2 Fragmentation Function and Jet Imbalance 

The key to measuring the fragmentation function of the jet of particles from 
a hard-scattered parton is to know the energy of the original parton which 
fragments, as pioneered at LEP [71] , Thus, in p-p collisions, a measurement 
of the direct-7 — h correlation from g + q — >> 7 + q, where the h represents 
charged hadrons opposite in azimuth to the direct-7, is (apart from the low 
rate) excellent for this purpose since both the transverse momentum and 
identity of the jet (8/1 ?/-quark, maybe 8/2 if the g + q 7 + g channel is 
included) are known to high precision. Two particle correlations are ana- 
lyzed in terms of the two variables [75]: p ou t = prsin(A0), the out-of-plane 
transverse momentum of a track; and where: 



z trig ~ PTt/PT]et is the fragmentation variable of the trigger jet, and z is 
the fragmentation variable of the away jet. Note that xe would equal the 
fragmenation fraction z of the away jet, for z t rig — 1, if the trigger and 
away jets balanced transverse momentum. The beauty of direct-7 for this 
purpose is that z tr ig = 1. 

Following the approach of Borghini and Wiedemann [72] who predicted 
the medium modification of fragmentation functions in the hump-backed or 
£ = ln(l/z) representation, PHENIX measured xe distributions in p-p col- 
lisions [73] and converted them to the £ = — In xe representation (Fig. |22^ t) 
which are in quite excellent agreement with the dominant ix-quark fragmen- 
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Figure 22: a) (right) £ = — In xe distributions for PHENIX isolated direct- 
7 p-p data [73] for all pr t ranges combined, compared to e + e~ collisions 
at = 14 and 44 GeV. b) (left) £ = -In distribution in Au+Au 

(0-20%) eu 

tation functions measured in e + e~ collisions at y/s/2 = 7 and 22 GeV, 
which cover a comparable range in jet energy. PHENIX preliminary mea- 
surements [H] of the £ = —Iuxe distribution in central (0-20%) Au+Au 
collisions, which suggest a modification consistent with Ref. [72] are shown 
in Fig. [22) 3. A firm conclusion awaits final results with improved statistics. 
Fragmentation functions from full jet reconstruction in A+A collisions are 
not yet available at RHIC. 

One of the important lessons learned at RHIC [75] about fragmentation 
functions is that the away-side xe distribution of particles opposite to a 
trigger particle (e.g. a 7r°), which is itself the fragment of a jet, does not 
measure the fragmentation function, but, instead, measures the ratio of pT a 
of the away-parton to pT t of the trigger-parton and depends only on the 
same power n as the invariant single particle spectrum: 



dP 



dxi 



N(n 



1 



x h U + ff) n 



(5) 



This equation gives a simple relationship between the ratio, xe ~ PT a /PT t = 
zt, of the transverse momenta of the away-side particle to the trigger parti- 
cle, and the ratio of the transverse momenta of the away-jet to the trigger-jet, 
Xh — PTa/ PT t - PHENIX measurements [76J of the xe distributions of 7r°-h 
correlations in p-p and Au+Au collisions at ^snn = 200 GeV were fit to 
Eq. [5] (Fig. [23^ i) [77J . The steeper distribution in Au+Au shows that the 
away parton has lost energy. The results for the fitted parameters are shown 
on the figure. In general the values of x 9 ^ do not equal 1 but range between 
0.8 < Xy- < 1.0 due to fey smearing and the range of xe covered. In order to 
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Figure 23: (left) a) xe distributions [77] from p-p (circles) and AuAu 0- 
20% centrality (squares) for pT t = 7 — 9 GeV/c, together with fits to Eq. [5] 
(solid lines) with parameters indicated. The ratios of the fitted parame- 
ters for AuAu/pp are also given, b) (right) Fractional jet imbalance [77] . 
1 - xfi A /x p h p , for RHIC and CMS data. 

take account of the imbalance {x p ^ < 1) observed in the p-p data, the ratio 
x AA /x vp is taken as the measure of the energy of the away jet relative to the 
trigger jet in A+A compared to p-p collisions. The fractional jet imbalance 
was also measured directly with reconstructed di-jets by the CMS collabo- 
ration at the LHC in Pb+Pb central collisions at y^NN — 2.76 TeV [66] . 
They also showed a large effect in p-p collisions. I calculated x^ from their 
p-p and Pb+Pb results [77] and compared these LHC values of 1 — x AA jx vp 
to those from PHENIX (Fig. (23)3). The large difference in fractional jet 
imbalance between RHIC and LHC cm. energies could be due to the dif- 
ference in jet p Tt between RHIC (~ 20 GeV/c) and LHC (~ 200 GeV/c), 
or different sensitivity of the direct and indirect methods, or the difference 
in n for the different y/s, or to a difference in the properties of the medium. 
Future measurements will need to sort out these issues by extending both 
the RHIC and LHC measurements to overlapping regions of pt- 

8 Anisotropic flow (v?) of direct-7 

Although direct-7 production [65] is the most beautiful QCD subprocess, it 
has a very serious problem: an overwhelming background of photons from 
higher tt° 7+7 and 77 — >> 7+7 decays makes it a very difficult experiment. 
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One must measure all the background sources: tt°, 77, ... , and calculate their 
contributions to the inclusive 7-ray spectrum. In principle, the background 
can be calculated whatever the pt distribution of the tt° and 77. However 
nature has been kind in that the invariant cross section for hard-scattering 
is a power law, da/pxdpT oc with n = 8.1 ± 0.05 for tt° with pt > 3 



GeV/c (Fig. |21Jl) [68]; also for p T > 3 GeV/c, ///tt = 0.48 ± 0.03 is a 
constant. This implies that for tt° 7 + 7, the spectrum of decay photons 
has the same power as the parent tv° so that the ratio at any pt is a constant: 

(6) 



7T 



T o n ■ 



The resulting background inclusive 7 spectrum from tv° and 77 decays at 
^/s^ = 200 GeV is: 

7back g roundA° « (1 + 0.48 x 0.39) x 2/7.1 = 1.19 x 2/7.1 = 0.334 (7) 

where 0.39 is the branching ratio for 77 — >► 7 + 7. In PHENIX we plot what 
we call the double ratio: 

-R7 ~ (Tinclusive/^ )/(Tbackground/^ ) ~ Tinclusive/Tbackground 

where it is important to see the calculated 7backgroimd/^° ratio, which usually 
comes from some opaque Monte Carlo program, to understand whether it 
makes sense according to Eqs. [6]and[7| 



Fig. 24=i [79J shows i? 7 for real photons, measured in an EM calorimeter, 
and virtual photons, which are e + e _ pairs from internal conversion of the 
direct-7, with 0.12 < m ee < 0.30 GeV/c 2 where there is no background from 
7T° Dalitz decay. This reduces the background by a factor of > 1.19/0.19 « 
6 [27], and allows the precision of i? 7 to be greatly improved as shown. 
Then, using the precise virtual photon i? 7 with the much higher statistics 
inclusive real-7 data, one can derive V2 for direct-7 from the measured V2 of 
inclusive real-7 compared to the measured V2 of 7's from tt° and rj decay. 

The result (Fig. [24)3) [79] is that the V2 of direct-7 is large in the range 
1 < Pt ^ 3 GeV/c but drops to zero for pt > 5 GeV/c where the photons are 
produced from initial hard-scattering and do not interact with the medium 
so that they do not flow. Since thermal radiation is produced in the flowing 
medium, the observed large V2 in the range 1 < pt < 3 GeV/c, where the 
direct-7 pt spectrum is exponential (Fig. pZlfc), confirms that these 7 are 
thermal radiation from the medium. What is very surprising is that the V2 
of the thermal photons is so large, the same or slightly greater than that of 

7T°'S. 
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Figure 24: a) (left) i? 7 vs pt [79J for virtual photons (solid circles) [69] 
and real photons (open squares) [78] for minimum bias Au+Au collisions at 
^/snn — 200 GeV. b) (right) Ratio of direct-7 V2 to tv° V2 for two reaction 
plane detectors [79] . 



9 A charming surprise 

This past year, the ALICE experiment at LHC [80] confirmed, with recon- 
structed charm mesons, the suppression of heavy quarks comparable to that 
of TV (from light quarks) for pr^4 GeV/c as previously observed at RHIC 



using direct-single-e^ from heavy quark (c, b) decay (Fig. 25 a) [81]. Also 
seen at RHIC is that heavy quarks exhibit collective flow (i^) (Fig. [25)3), 
another indication of a very strong interaction with the medium. The fact 
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Figure 25: PHENIX |H]: a) (left) R A a (central Au+Au) b) (right) v 2 
(minimum bias Au+Au) as a function of pt for direct-e 1 * 1 at ^s^n — 200 
GeV. 
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that heavy quarks are suppressed the same as light quarks strongly disfavors 
the QCD energy-loss explanation of jet-quenching because, naively, heavy 
quarks should radiate much less than light quarks and gluons in the medium; 
but opens up a whole range of new possibilities including string theory [82] . 

10 Zichichi to the rescue? 

In September 2007, I read an article by Nino, "Yukawa's gold mine" in the 
CERN Courier taken from his talk at the 2007 International Nuclear Physics 
meeting in Tokyo, Japan, in which he proposed: "We know that confinement 
produces masses of the order of a giga-electron-volt. Therefore, according 
to our present understanding, the QCD colourless condition cannot explain 
the heavy quark mass. However, since the origin of the quark masses is still 
not known, it cannot be excluded that in a QCD coloured world, the six 
quarks are all nearly massless and that the colourless condition is 'flavour' 
dependent." 

Nino's idea really excited me even though, or perhaps because, it ap- 
peared to overturn two of the major tenets of the Standard Model since it 
seemed to imply that: QCD isn't flavor blind; the masses of quarks aren't 
given by the Higgs mechanism. Massless b and c quarks in a color-charged 
medium would be the simplest way to explain the apparent equality of gluon, 
light quark and heavy quark suppression indicated by the equality of Raa 
for 7T° and direct single-e^ in regions where both c and b quarks dominate. 
Furthermore RHIC and LHC-Ions are the only place in the Universe to test 
this idea. 

Nino's idea seems much more reasonable to me than the string theory 
explanations of heavy-quark suppression (especially since they can't explain 
light-quark suppression). Nevertheless, just to be safe, I asked some distin- 
guished theorists what they thought, among others [27], Steve Weinberg. 
He said that he and Lenny Susskind had a model Technicolor (or Hyper- 
color) that worked well in the vector boson sector but didn't give mass to 
the fermions. 

Nino proposed to test his idea by shooting a proton beam through a QGP 
formed in a Pb+Pb collision at the LHC and seeing the proton 'dissolved' 
by the QGP. My idea is to use the new PHENIX VTX detector, installed 
in 2011, to map out, on an event-by-event basis, the di-hadron correlations 
from identified b — b di-jets and identified c — c di-jets, which do not originate 
from the vertex, and light quark and gluon di-jets, which originate from the 
vertex and can be measured with 7r°-hadron correlations. A steepening of 
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the slope of the xe distribution of heavy-quark correlations as in Fig. [23^ t 
will confirm in detail (or falsify) whether the different flavors of quarks 
behave as if they have the same energy loss (hence mass) in a color-charged 
medium. If Nino's proposed effect is true, that the masses of fermions are 
not given by the Higgs particle and all quarks are nearly massless in the 
QGP, and we can confirm the effect at RHIC or LHC-Ions, this would be a 
case where we Relativistic Heavy Ion Physicists may have something unique 
to contribute at the most fundamental level to the Standard Model, which 
would constitute a "transformational discovery." Of course the LHC could 
falsify this idea by finding the Higgs decay to b — b at the expected rate in 
p-p collisions. Clearly, there are exciting years ahead of us! 
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